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ABSTRACT 

I We consider the relic neutrino background produced by Population III stars cou- 

pled with a normal mode of star formation at lower redshift. The computation 
is performed in the framework of hierarchical structure formation and is based 



o 



' on cosmic star formation histories constrained to reproduce the observed star 

formation rate at redshift z < 6, the observed chemical abundances in damped 
Lyman alpha absorbers and in the intergalactic medium, and to allow for an 
0^ I early reionization of the Universe at z ^ 10 — 20. We consider both a burst and 

non-burst model for Population III star formation. We find that although the 
, high redshift burst of Population III stars does lead to an appreciable flux of 

\^ • neutrinos at relatively low energy {E^ « 1 MeV), the observable neutrino flux is 

I dominated by the normal mode of star formation. We also find that predicted 

CN| ' fluxes are at the present level of the SuperK limit. As a consequence, the su- 

pernova relic neutrino background has a direct impact on models of chemical 
' evolution and/or supernova dynamics. 

o 

1. Introduction 

o 

+i . Among the resuhs contained in the first year data obtained by WMAP^'-' was 

. large optical depth implying that the Universe became reionized at high redshift in 

the range, 11 < z < 30 at 95% CL. This can be explained by a generation of massive 
stars formed at high redshift!^. Such a scenario would have notable consequences 



X 

^ I for cosmic chemical evolution, in particular, the metal enrichment of the interstellar 

medium (ISM) and intergalactic medium (IGM) '^ ' '^ ' ^ . In addition, the cosmic star 
formation rate (SFR) at redshifts 2; < 6 is observed to be over an order of magnitude 
larger than the current SFR and is seen to have peaked at redshift z ~ 3 

A direct consequence of this is the predicted enhancement in the rate of core 
collapse supernovae. In addition to the sharp spike of supernovae at very high redshift 
due to the explosions of stars responsible for the early epoch of reionization, the 
enhanced SFR of the normal mode of star formation at redshifts z < 6 leads to a 
supernova rate which is approximately a factor of 30 times the current rate, and a 
factor of 5 times the observed rate at z ~ 0.7 

Another consequence of an enhanced SFR and SN rate is the resultan t neutrino 
background spectrum produced by the accumulated core collapse 

supernovadHMmmniini. 

^Summary of talk given at the Ilird International Workshop on: NO-VE "Neutrino Oscillations 
in Venice" , Venice Italy, February 2006 



Here, we incorporate fully developed chemical evolution models which trace the his- 
tory of pre-galactic structures and are based on a ACDM cosmology and include a 
Press-Schechter modei^ of hierarchical structure formation. This will allow us to 
make detailed estimates of the predicted fluxes of neutrinos which can be detected in 
large underground detectors. 

2. Chemical Evolution 

The basic elements of any chemical evolution model include: the star formation 
rate, ip{t), which is observed to be greatly enhanced at high redshift; the initial mass 
function, (p{m), assumed to have a near Salpeter slope of 1.3; element yields to allow 
us to compare with the chemical abundances in DLAs, the IGM and metal-poor stars. 

The model employed is an open-box type model, which includes infall due to the 
accretion of baryons in the hierarchal scenario, and outflows of gas due to supernova 
winds. The cosmic star formation histories we consider here have been adopted from 
detailed chemical evolution models^ These models are bimodal and are described 
by a birthrate function of the form 



where 0i(2) is the IMF of the normal (massive) component of star formation, and ipi(2) 
is the respective star formation rate. Z is the metallicity. The normal component 
contains stars with masses between 0.1 and 100 Mq and is primarily constrained by 
observations at low redshift (2; < 6). The massive component operates at high redshift 
and is cut off once the metallicity reaches a critical value taken to be 10~^. Both 
components can contribute to the chemical enrichment of galaxy forming structures 
and the IGM. We examine three different models of the massive mode as described 
below. 

We consider several minimal masses for minihalos for star formation to occur: 
10^, 10^, 10^, 10^, and 10^^ M0. Star formation is assumed to begin when the baryon 
fraction in the structures reaches /{, = 0.01. This criteria fixes the initial redshift for 
star formation to begin. For example, for = 10'^ Mq, star formation begins at 

z = 16. For larger (smaller) masses, star formation begins later (earlier). Later we 
will relax this assumption. 

3. The normal mode 

The normal mode of star formation is referred to as Model and provides a 
standard star formation history, with stellar masses in the range 0.1 M© < m < 100 
Mq. The SFR for Model is of an exponential form 



B{m, t, Z) = 0i(m)^i(t) + ^2{m)^2{Z) 



(1) 



■t/n 



(2) 
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Figure 1 : The cosmic star formation rate for the normal mode of star formation (Model 0) . Plotted 
are the results of the best fit models for each choice of Afmin- 



which corresponds to a SFR dominated by elhptical galaxies. Best fits for ui and ti are 
given in Table 1. In thes e mo dels the outfiow is non-zero and the details for computing 
the outfiow are given in 1^1^. The overall efficiency of outfiow is parameterized by e 
whose value is also given in Table 1. This model alone is inadequate for high redshift 
reionization. 
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Table 1: The model parameters for the normal mode of star formation (Model 0). Column 1 indicates 
the input value of the minimum mass for star forming structures. Column 2 is derived from column 
1, having assumed that = 1% when star formation begins. In columns 3, 4, and 5, parameter 
values for the efficiency of outflow and the SFR are given. The slope of the IMF is x = 1.3 for all 
models. 

The results of the best fit models are plotted in Figure [T] for each value M^i^i as 
indicated. Because we have fixed the initial baryon fraction in structures, each value 
of Mmin corresponds to a different initial redshift. As one can see, each of the models 
gives a satisfactory fit to the global SFR, save perhaps the case with M-^^^ = 10^^ 
Mq. The data shown (taken from^) has already been corrected for extinction. 

Some of the consequences of the SFRs shown in Figure Q are shown in Figure El 
where the type II supernova rate and metal enrichment are shown. The observed rates 
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Figure 2: The type II supernova rate and global metallicity for Model 0. Plotted are the results of 
the best fit models for each choice of Mmin- 

of type II supernovae up to 2; ~ 0.7 is taken froml^. The predicted type II supernova 
rates are also an excellent fit to the existing data at low redshift. The predicted 
rates are substantially higher at higher redshift and lower Mmin- Also shown is the 
evolution of the metallicity in both the ISM and IGM (dashed curves). While the ISM 
metallicity rises very quickly initially, evolution is more gradual in the IGM where 
it is controlled by the outflow. The metallicity of 100 Damped Lyman-a systems as 
measured by is also plotted. Note that only the normal mode is included here, 
and contributions from Population III stars will only serve to increase the metallicity 
in both the ISM and IGM. 

4. The massive mode 

We consider three different models, labeled Models 1, 2a, and 2b to describe the 
massive mode. They are distinguished by their respective stellar mass ranges. In 
Model 1, the IMF is defined for stars with masses in the range 40 Mq < m < 100 
Mq. All of these stars die in core collapse supernovae leaving a black hole remnant 
with the mass of the progenitor's helium core. They all contribute to the chemical 
enrichment of the ISM and IGM. This period of star formation is brief and is described 



by a SFR of the form 

= ^^2M,tructe-^^°^/^-'' (3) 

where Mgtruct is the mass of the primitive star forming structure and Zrr t = 10~^ 
Zq is the critical metalhcity at which Population III star formation ends^. Other 
parameters describing these models are given in Table 2 for the case of M^[^ = 10^ 
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Table 2: Parameter values for the massive starburst Models 1, 2a and 2b. Column 1 indicates the 
model number and column 2 the input value of the minimum mass for star forming structures. In 
columns 3 and 4, we show the adopted outflow efflciency and the coefficient of the massive mode 
SFR. 

When a massive mode is added to the normal mode described by Model 0, the 
outflow efficiency must be adjusted so as to avoid the overproduction of metals in the 
IGM. However, there is a degeneracy in the massive mode parameters e and z/2. In the 
model labeled 1, the massive mode contributes roughly 50% of the IGM metallicity 
at a redshift z = 2.5. By increasing z/2 and decreasing e, this contribution can be 
increased to 90%, at the same time increasing the ionization capacity of the model. 
This case is labeled le. 

As seen in Figure El star formation in Model 1 begins at a very high rate and 
falls precipitously as metals are injected into the ISM. Below we will also consider a 
variant of Model 1 in which the massive mode begins much earlier but its duration 
is extended in redshift. As in the case of Model 0, the onset of star formation is 
determined by Mmin and the initial value for the baryon fraction in structures (fixed 
to be 1%). As one can see in the lower panel, the metallicity in the ISM reaches 
values far in excess of Z^it due to the finite lifetime of the massive stars relative to 
the speed at which the metallicity is attained. Notice also that once the metallicity 
from Pop III stars is produced, the ISM metallicity later decreases as a result of the 
accretion of metal-free g the structures grow. 

Among the chief motivating factors in developing a model of cosmic chemical 
evolution is the early reionization of the Universe. Several studies suggest that an 
early burst of star formation, as in Model 1, is sufficient? -'^ ^ In Figure EJ we 
show the number of ionizing photons per baryon produced for each of our choices of 
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Figure 3: The cosmic star formation rate and global metallicity for the normal mode with an added 
burst of massive stars (Model 1). 



Miain- The procedure for computing this stellar ionizing flux is explained inl^. It is 
important to remember that only a fraction /esc of these UV photons will escape the 
structures and therefore be available to ionize the IGM. The effective value of /esc is 
poorly known but could vary from about 1 to 30%. The minimum number of photons 
required for complete reionization is also plotted for three possible dumpiness factors. 
The ionizing potential clearly increases with Mmin and decreasing redshift. The ratio 
of this minimum number (dashed line) to that for the stellar ionizing photons (solid 
line) gives the minimum fraction /esc necessary to fully reionize the IGM. 

Model 2a is described by very massive stars which become pair instability super- 
novae. The IMF is defined for 140 Mq < m < 260 Mq and the SFR for this model is 
the same as that for Model 1, but must be reduced by a factor of 8 due to constraints 
on metal abundances in the ISM. The most massive stars are considered in Model 2b 
and fall in the range 270 Mq < m < 500 Mq, with the SFR as in Model 1. These 
stars entirely collapse into black holes and do not contribute to the chemical enrich- 
ment of either the ISM or IGM. A comparison of the star formation rates, ionization 
potentials, and global metallicity in models 1, 2a, and 2b is shown in Figure El 
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Figure 4: The cumulative number of ionizing photons emitted by stars is plotted as a function of 
redshift for all models of Figure |3| The minimum number of photons per baryon necessary to fully 
reionize the IGM is plotted as a function of redshift by thin dashed lines for three different values 
of the dumpiness factor Ch ii — 10, 30 and 100. 

Rather than the rapid burst models for Population III stars discussed above, it 
is useful to consider the consequences of an extended period of Population III star 
formation as described in 1^1. The non-burst models begin at high redshift, z ^ 30, 
independent of the the minimum halo mass. As such, they do not depend on any 
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Figure 5: The star formation rate, reionization potential, and global metallicity for Models 1, 2a, 
and 2b with A^min = 10^ Mq. 

choice of the initial baryon fraction. In this case, we adopt a SFR given by 

^2 = i/2MisMe-^^°^/^"'' (4) 

Notice the replacement of the structure mass with the mass of gas in the structure 
ISM. 

Given the new SFR, it is necessary to refit the observed data and once again 
maximize z/2 for ionization without the overproduction of metals. New parameter 



values of outflow and the SFR are given in Table 01 for non-burst Models 1, 2a, and 
2b, now labeled In, 2an, and 2bn respectively. 
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Table 3: As in Table |21 parameter values for nonburst Models In, 2an and 2bn. 



In Figure El we show the SFR, ionization potential, and metallicity for such a non- 
burst model compared with the analogous result for Model le, restricting attention 
to Mjnin = lO'^ Mq. In this case, massive star formation begins at z = 30. The 
dotted curve in the top panel of Figure IHl is the massive component of Model In. 
While the massive mode SFR in Model In is significantly lower than the peak burst 
of le, its duration is far longer and continues down to redshifts < 3, although the 
total star formation rate becomes dominated by the normal mode when z < 10. The 
clear difference between Model In and the rapid burst Model le is the production of 
ionizing photons at high redshift as shown in the middle panel. The burst is capable 
of ionizing the IGM at z = 16 for all three choices of Ch h. In contrast, the metallicity 
shown in the lower panel is relatively unaffected by the rapid burst. 

Finally, given an IMF and SFR, it is straightforward to compute the rate of core 
collapse supernovae, 

SNR= / dm (l){m)ij {t - T{m)) , (5) 

where mjnin(^) is the minimum mass with lifetime, r, less than t. 



5. Neutrino Production 



In all of the models we will consider, star formation begins at high redshift, domi- 
nated initially by massive stars which may explode as core collapse or pair-instability 
supernovae and provide for the reionization of the IGM. Each explosion, regardless 
of type, releases most of the star's gravitational energy in the form of neutrinos with 
a specified energy spectrum and flux. Given a chemical evolutionary model, the 
integrated contribution of SN to the neutrino background can be computed. 

The expected differential flux of neutrinos at Earth with energy E can be expressed 

as 
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Figure 6: The star formation rate, ionizing photons per baryon and global metallicity for Model In 
and a rapid burst Model le with M^in ~ 10'' Mq. The dotted line in the upper panel shows the 
SFR of the massive mode of Model In. 



where N^^{m) is the total number of neutrinos of a given species, a, emitted in the 
core collapse of a star of mass m, and represents the neutrino energy spectra at 
emission with energy E' = E{1 + z). The integration limits M^m and M^ax are the 
minimum and maximum masses in each model for which supernovae occur, and Zi is 
the initial redshift at which star formation begins. 

When a star undergoes core collapse, the mass of the remnant is determined by 
the mass of the progenitor. We assume that all stars of mass m > 8 Mq will die as 
supernovae. For stars of mass 8 Mq < m < 30 Mq, the remnant after core collapse 
will be a neutron star of m ~ 1.5 Mq. More massive stars fall into two categories; 
black holes and pair instability supernovae. Pair instability supernovae are thought 
to occur for stars with 140 Mq < m < 260 Mq, in which case the explosion leaves no 
remnant. All other stars collapse to form black holes. Stars with 30 Mq < m < 100 
Mq become black holes with mass approximately that of the star's helium core before 
coUaps^^ . We take the mass of the Helium core to be 

13 

Mne = ^-{m- 20 Mq) (7) 



for a star with main sequence mass m We assume that stars with m > 260 Mq 
collapse entirely to black holes. 



The energy emitted in each core collapse, Ecc corresponds to the change in gravi- 
tational energy, 99% of which is emitted as neutrinos^. In the cases where collapse 



irrespc 

M. 

results in a neutron star, Ecc = 5 x 10^^ ergs. For stars that collapse to black holes, 
Ecc is proportional to the mass of the black hole. For masses less than 100 Mq, we 
take Ecc = O.SMhc ■ Pair instability supernovae experience a much more powerful 
explosion than core collapse supernovae, however few neutrinos are emitted and with 
very low energies such that they would not be observeci^. In this case, we assume 
the energy emitted in neutrinos is the same as that for ordinary core collapse super- 
novae yielding a neutron star, but the average neutrino energy is {EpJ = 1.2 MeV'^. 
For the most massive stars considered, Ecc = 0.3m. Although several studies find a 
distinct hierarchy in the partitioning of neutrino luminosity among the species dur- 
ing the different luminosity phases of cor e collaps e, equipartition of the total energy 
emitted by the star is generally acceptec0^^l^ . For a comparison of luminosity 
hierarchies found in recent simulations, see Keil et al.^^. 

We assume that each electron neutrino carries an average energy (Ei,^) = 13.3 MeV, 
which is consistent with simulations. The charged current reactions that prevent neu- 
trinos from emerging from the star are uqu — >■ pe~ and z/eP — >■ ne"*". The different 
trapping reactions result in different neutrinosphere radii, and therefore different aver- 
age energies for z/g and Ue- We assume {EpJ = 15.3 MeV. The other species, denoted 
u^, undergo only neutral current interactions. The mechanism that governs their 
average temperature at emission is more complicated, but the generally accepted hi- 
erarchy is {E^^) < {EpJ < {E„J. We have taken {E„J = 20 MeV. The total number 



of z/q, emitted by a star during core collapse is given by 



AT,, 



Ecc 



The neutrino spectra at emission can be described by a normalized Fermi-Dirac 
distribution, 



dPg 

dE' 



E' 



where = 180Cs{E,^J /Ttc'^ is the effective neutrino temperature taken to be inde- 
pendent of the mass of the star. We assume a flat ACDM cosmology with 
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dz 



9.78/i-iGyr 



(10) 



where Vt^ = 0.73, = 0.27, and h = 0.71^. 



6. Resulting Neutrino Fluxes 

Having described the chemical evolution models and the calculation of the neutrino 
spectrum, we now show the results for each of the models considerecA^ . Model 
fluxes are plotted in Figure [7| Results are shown for several choices of minimum halo 
masses, Mmin. As we will see below these fluxes are large enough to be probed by 
current detectors. As M^[^ is increased, star formation occurs at later redshift and 
as a result, the peak of the neutrino flux is shifted slightly to higher energy. 




2 4 6 8 

(MeV) 



Figure 7: Fluxes from Model for five choices of M^in = 10^ (sofid), lO'^ (thick), 10** (dotted), 10^ 
(dashed) and lO" (dot-dashed) Mq. 

The massive modes of Model 1 and le fluxes are plotted in Figure|H]for the specific 
choice of Mmin = 10'' Mq which is the preferred case in*^. As seen in FigureEl massive 



stars associated with Population III turn on at a redshift of approximately 16, but the 
duration of the burst is relatively brief. As a result, the peak of the flux distribution 
is at relatively low energy. More importantly, because of the brevity of the burst, 
the entire neutrino spectrum is redshifted down, in contrast to the Model spectrum 
which extends to higher energy due to stars produced at lower redshifts. As expected, 
the more extreme model, le, has a peak flux which is about 5 times that found for 
Model 1. This is directly related to the increased SFR in Model le as characterized 
by the increase in z/2. 
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Figure 8: Fluxes from the massive modes of Model 1 (black) and le (grey). 

Similarly, we show in Figure IHl the resulting flux from very massive Population III 
stars corresponding to Models 2a, 2ae and 2b''. As before, the fluxes from Models 
2a are relatively small and peak at very low energy as seen in the insert to the 
figure. In Figure ITUl we show the total fluxes in Models 1, 2a, and 2b with Mmin = 
10^ Mq. As one expects, the low energy spectrum is dominated by neutrinos produced 
in the massive mode, whereas the spectrum at higher energies {E,^ > 3 MeV), is 
indistinguishable between the models and dominated by the normal mode. 

Finally, we show the results for the non-burst Models In, 2an, and 2bn as com- 
pared with the burst models described above. In Figure ITT| we show the fluxes from 
the massive modes of models In and 2an compared with models 1 and 2a. As can 
be expected, the fluxes extend to higher frequencies due to the longer duration of 
the stellar population. The total fluxes including the normal mode are shown in Fig- 
ure El Similarly, in Figure El we compare the fluxes of Models 2b and 2bn. Because 
the fluxes are so large, they dominate the normal mode flux (over the energy range 
shown) and we show only the total flux for this case. 

Since stars associated with Model 2b do not contribute to element enrichment, there is no Model 

2be. 




7. Detectability 

The solar neutrino flux at Earth is larger than the expected flux from SRN's by 
several orders of magnitude for E^, < 19 MeV However, neutrinos, rather than 
antineutrinos, are produced in the thermonuclear reactions in the sun and have a 
smaller cross section for detection by about two orders of magnitude. This, with the 
directional information from recoil electrons in the detector, allows this background 
to be excluded at SuperK, KamLAND, and SNO. But decays of spalled nuclei from 
cosmic ray muons constitute an unavoidable background in this range. The current 
upper limit for the flux of SRN's at SuperK is 1.2cm-^s-^ for > ig.SMeV^Sl. in 
addition, it has been pointed out recently that if the background analysis from 
SK is coupled with the sensitivity to electron neutrinos at SNO it will be possible to 




Figure 11: Fluxes from the massive modes of Models In and 2an (grey) and 1 and 2a (black). 
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Figure 12: Total fluxes from Models In and 2an (grey) and 1 and 2a (black). 

reduce the upper limit on the flux of electron neutrinos. SNO should be sensitive to 
a flux of 6 cm-^s"^ in the range 22.5 MeV< E^^ < 32.5 MeV. 
In Figure we show the observable flux 

F(^thresh) = / ^ 

as a function of detector threshold energy. While the fluxes are quite appreciable at 
low threshold energies, they in fact remain relatively high at larger energies due to the 
large SFR associated with Model 0. Indeed, in all of our models, our predicted flux 
above 19.3 MeV already exceeds the current bound of 1.2 cm~^s~^ from SuperK^^. 
The detailed flux predictions are given in Table |3J where we show the detectable 
flux for the viable energy windows at SK and SNO for both the burst and non-burst 
models. Although SRN's will likely not be seen at SNO given these flux levels, in 
many of our models the SK bound is saturated by the expected flux, indicating that 
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Figure 13: Total fluxes from Models 2bn (grey) and 2b (black). 

SRN's may be observed in the near future. Models labeled osc, include the effects of 
neutrino oscillations. 




5 10 15 20 25 30 

(MeV) 



Figure 14: Detectable fluxes from Model 1 with (dashed) and without (solid) oscillations as a function 
of neutrino energy threshold. 

Despite the large fluxes displayed in Table E] relative to the SK limilP^, one can 
not conclude that the models considered have already been excluded by experiment. 
There are of course many uncertainties built into our chemical evolution models as 
well as uncertainties in the adopted neutrino physics. These fluxes are very sensitive 
to our assumed average neutrino energy. Recall our adopted value for Ep^ is 15.3 MeV. 
In Fig. El we show the sensitivity of the flux above 19.3 MeV to the average neutrino 
energy. In order to satisfy the SK limit of 1.2 cm~^s~^, we would have to lower {EpJ 
to 13.3 MeV. This is fully consistent with the range of neutrino energies in supernova 
models . 
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Table 4: Predicted fluxes in cm ^ in the models considered here. Results are given for electron 
antineutrinos with energies E;^^ > 19.3MeV for SK and for electron neutrinos with 22.5MeV< E^^ < 
32.5MeV for SNO. 



8. Summary 

We have considered several scenarios for star formation which reproduce the ob- 
served chemical abundances and SFR for z < 6 and reionize the universe at high 
redshift. Each model of star formation here consists of a normal mode coupled to a 
Population III mode of massive star formation at high redshift. We considered both 
a burst and non-burst model for Population III star formation. 

In the burst model, because the massive mode of star formation is so brief and 
takes place at high redshift, the corresponding electron anti-neutrino fluxes peak at 

<^ 1 MeV. Even in the non-burst model, fluxes are become small at energies above 
~ 4 MeV. Thus despite the large fluxes produced by the massive mode, these low- 
energy neutrinos will be difficult to detect. In contrast, the normal mode of star 
formation, which dominates the flux at observable energies, is peaked at a somewhat 
higher energy and has a broad spectrum due to the production of stars at lower 
redshift. Our calculated fluxes of SRN's from core collapse, however, saturate the 
SK bound of 1.2cm^^s^^ for E^, > 19.3 MeV in all models. Although there are 
uncertainties in the neutrino physics, such as the average energies at emission, the 
prospects for observation in the near future are good. 
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